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Abstract 

We evaluate the Casimir energy and force for a massive fermionic field in the geometry 
of two parallel plates on background of Minkowski spacetime with an arbitrary number of 
toroidally compactified spatial dimensions. The bag boundary conditions are imposed on 
the plates and periodicity conditions with arbitrary phases are considered along the compact 
dimensions. The Casimir energy is decomposed into purely topological, single plate and 
interaction parts. With independence of the lengths of the compact dimensions and the 
phases in the periodicity conditions, the interaction part of the Casimir energy is always 
negative. In order to obtain the resulting force, the contributions from both sides of the 
plates must be taken into account. Then, the forces coming from the topological parts of the 
vacuum energy cancel out and only the interaction term contributes to the Casimir force. 
Applications of the general formulae to Kaluza-Klein type models and carbon nanotubes are 
given. In particular, we show that for finite length metallic nanotubes the Casimir forces 
acting on the tube edges are always attractive, whereas for semiconducting-type ones they 
are attractive for small lengths of the nanotube and repulsive for large lengths. 



PACS numbers: 03.70.+k, ll.10.Kk, 61.46.Fg 



1 Introduction 

A key feature of most high energy theories of fundamental physics, including supergravity and 
superstring theories, is the presence of compact spatial dimensions. From an inflationary point of 
view universes with compact spatial dimensions, under certain conditions, should be considered 
a rule rather than an exception [TJ. The models of a compact universe with non-trivial topology 
may play an important role by providing proper initial conditions for inflation (for physical 
motivations of considering compact universes see also [2]). There has been a large activity to 
search for signatures of non-trivial topology by identifying ghost images of galaxies, clusters 
or quasars. Recent progress in observations of the cosmic microwave background provides an 
alternative way to observe the topology of the universe [3]. If the scale of periodicity is close to 
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the particle horizon scale then the changed appearance of the microwave background sky pattern 
offers a sensitive probe of the topology. An interesting application of the field theoretical models 
with compact dimensions recently appeared in nanophysics [3j . The long- wavelength description 
of the electronic states in graphene can be formulated in terms of the Dirac-like theory in 3- 
dimensional spacetime with the Fermi velocity playing the role of speed of light (see, e.g., Refs. 
|5j). Single- walled carbon nanotubes are generated by rolling up a graphene sheet to form a 
cylinder and the background spacetime for the corresponding Dirac-like theory has topology 
R 2 xS\ 

In quantum field theory the boundary conditions imposed on fields along compact dimensions 
change the spectrum of vacuum fluctuations. The resulting energies and stresses are known 
as topological Casimir effect (for the topological Casimir effect and its role in cosmology see 
[6]- [H] and references therein). In the Kaluza-Klein-type models this effect has been used as 
a stabilization mechanism for moduli fields which parametrize the size and the shape of the 
extra dimensions. The Casimir energy can also serve as a model of dark energy needed for the 
explanation of the present accelerated expansion of the universe (see |12j and references therein) . 
In addition to its fundamental interest the Casimir effect also plays an important role in the 
fabrication and operation of nano- and micro-scale mechanical systems (see, for instance, |13j). 

The effects of the toroidal compactification of spatial dimensions on the properties of quan- 
tum vacuum for various spin fields have been discussed by several authors (see, for instance, 
[6]-[ll], [141 115} [TB] and references therein). The combined effect of extra compactified dimen- 
sions and boundaries on the Casimir energy in the classical configuration of two parallel plates 
has been recently considered in [17] for a scalar field and in [18] for the electromagnetic field. 
The Casimir energy and forces in braneworld models have been evaluated in Refs. [19] by using 
both dimensional and zeta function regularization methods. Local Casimir densities in these 
models were considered in Refs. |20j . The Casimir effect in higher dimensional generalizations of 
the Randall-Sundrum model with compact internal spaces has been investigated in [21]. In the 
present paper, we investigate the Casimir effect for a massive fermionic field in the geometry of 
two parallel plates on background of spacetime with an arbitrary number of toroidally compact- 
ified spatial dimensions. We will assume generalized periodicity conditions along the compact 
dimensions with arbitrary phases and MIT bag boundary conditions on the plates. This problem 
in background of 4-dimensional Minkowski spacetime with trivial topology has been considered 
in [22] for a massless field and in [23] in the massive case (see also [6]). For arbitrary number 
of dimensions the corresponding results are generalized in Refs. [24} [25] for the massless and 
massive cases respectively. The Casimir problem for fermions coupled to a static background 
field in one spatial dimension is investigated in [26]. The interaction energy density and the 
force are computed in the limit that the background becomes concentrated at two points. The 
fermionic Casimir effect for parallel plates with imperfect bag boundary conditions modelled by 
5-like potentials is studied in [27] . 

This paper is organized as follows. In the next section, we specify the eigenfunctions and the 
eigenmodes for the Dirac equation in the region between the plates assuming the bag boundary 
conditions on them. In section [31 by using the Abel-Plana-type summation formula, we present 
the Casimir energy in the region between the plates as the sum of pure topological, single 
plate and interaction parts. In section 2] we consider the Casimir force acting on the plates. 
In section [5] we evaluate the Casimir energy and forces by making use of the generalized zeta 
function technique. An alternative representation of the single plate part of the Casimir energy 
is also given. The special case of topology R D ^ X x S 1 is discussed in section [6] In section [7] we 
give applications of general formulae to the Casimir effect for electrons in finite-length carbon 
nanotubes within the framework of 3-dimensional Dirac-like model. The main results of the 
paper are summarized in section [8] 
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2 Eigenfunctions and eigenmodes 

We consider a quantum fermionic field tp on background of {D + l)-dimensional flat spacetime 
with spatial topology R p+1 x (S" 1 ) 9 , p + q+ 1 = D. The corresponding line element has the form 

D 

ds 2 = dt 2 - ^(d^) 2 , (1) 
i=i 

where — oo < z l < oo, I = 1, . . . ,p + 1, and ^ z ^ L\ for / = p + 2, . . . , D. We assume that 
along the compact dimensions the field obeys boundary conditions 

i>(t, z p , z p+ \z q + L iei ) = e 2ma '^(t, zp, z,), (2) 

with constant phases ^ a\ < 1. In (J2j) , z p = (z 1 ,...,z p ) and z,j = {z p+2 , . . . , z D ) denote 
the coordinates along uncompactified and compactified dimensions respectively, e; is the unit 
vector along the direction of the coordinate z l , I = p + 2, . . . ,D. The periodicity conditions 
for untwisted and twisted fermionic fields are obtained from (|2|) as special cases with oti = 
and ai = 1/2 respectively. As we will see below, special cases a\ = 0,1/3,2/3 are realized in 
nanotubes. 

In this paper we are interested in the Casimir effect for the geometry of two parallel plates 
placed at z p+l = and z p+1 = a on which the field obeys the MIT bag boundary condition: 

(1 + i^n^) i/j = , z p+1 = 0, a, (3) 

where 7^ are the Dirac matrices and is the normal to the boundaries. In the (D + 1)- 
dimensional spacetime the Dirac matrices are N x N matrices with N = 2^ D+l ^ 2 \ where the 
square brackets mean the integer part of the enclosed expression. We will assume that these 
matrices are given in the chiral representation: 

with the relation 0^0^ + cr u a^ = 25 In the discussion below we consider the region between 

the plates, ^ z p+1 ^ a, where we have n M = — <5^ +1 for the plate at z p+1 = and n M = Sf^ 1 
for z p+1 = a. 

The dynamics of the field is governed by the Dirac equation 

ij^d^ip — mip = . (5) 

Assuming the time dependence in the form e i£Jt , the positive- and negative- frequency solutions 
to this equation can be presented as 

4+) = Aae-**( . . A 
y p p \ -ia+ -V(p/ (to + m) /' 

= v^(^' Vx/ ^ +m A (6) 

where a = (<ti, . . . ,<td), w = y^k 2 , + kp+i + k 2 + m 2 , and 

<p = e ik ir z ll L +e ik ^ lzP+1 +^_e-^+ lzP+1 ) , 

x = e - ik ir z ll ( x+e ^+i^ +1 +X-e - ^ +l2!P+1 ) , (7) 
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with k|| = (k p , k 9 ) and k p = (k\, . . . , k p ), k g = (fc p+2; • • • > The eigenvalues for the compo- 
nents of the wave vector along the compactified dimensions are determined from the periodicity 
conditions 



k <? = ( 27r ( re p+2 + a p+2 )/L p . 



+2, 



,2Tr(n D + a D )/L D ) 



n 



p+2, 



,n D = 0,±l,±2,.... (8) 



For the components along the uncompactified dimensions one has — oo < ki < oo, I = 1, . . . ,p. 

From the boundary condition on the plate at z p+1 = we find the following relations between 
the spinors in © 



m(to + m) + k? 



X- 



p+i 



kp+iVp+io-f ■ k|| 



(m 

m(u) + m) + k 2 ^ — kn^rr . , , cm 



ikp+i) (uj + m) 



■<P- 



-x+, 



(9) 



where cry 
with 



(m + i/c p+ i)(u; + m) 
,<7d). We will assume that they are normalized in accordance 



xXx+ = i. 

As a set of independent spinors we will take (p- = and xh 



(10) 



( CT )' , where 
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1, . . . , iV/2, are one-column matrices having iV/2 rows with the elements uj, = 5i a , and w^' = 
iw( a \ Now the set of quantum numbers specifying the eigenfunctions © is (3 = (k, a). From 

a it follows that the eigenvalues of k p+ \ are roots of the 



the boundary condition at z p+1 
transcendental equation 



ma sin(k p+ ia) / (k p+ ia) + cos(k p+ ia) 



0. 
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All these roots are real. We will denote the positive solutions of Eq. (jilt) by A ra = k p+ ia, 
n = 1,2,.... For a massless field we have A n = ir(n — 1/2). Note that the equation (jlip 
determining the eigenvalues for /c p +i does not contain the parameters of the compact subspace 
and is the same as in the corresponding problem on the topologically trivial Minkowski spacetime 
(see [6]). 

The normalization coefficient A p in ([6]) is determined from the orthonormalization condition 



dza 



' dz p+1 ^ ±)+ # 



y /3' 



(12) 



where 5ppi is understood as the Dirac delta function for continuous indices and the Kronecker 
delta for discrete ones. Substituting the eigenfunctions ([6]) into this condition one finds 



A, 



m 



4(2ir)PuaV q 



1 



sin(2A;, 



p+1' 



2k p j r \a 



(13) 



where V q = L p+2 ■ • • Ad is the volume of the compact subspace. 



3 Casimir energy 

For the spatial topology R p+1 x (S* 1 ) 9 the vacuum energy (per unit volume along the directions 
z z p ) in the region between the plates is given by the following mode-sum: 

J V ; n,eZH!=l 
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where n g = (n p+1 , ...,«£>) and 



u, 2 



+ k 2 nq + X 2 Ja 2 + m 2 , k 2 nq = £ [2^(n, + a/ )/L,] 2 . (15) 

«=p+2 



Of course, the expression on the right hand-side of Eq. (|14p is divergent. We will assume that 
some cutoff function is present, without writing it explicitly. For the further evaluation of the 
Casimir energy we apply to the sum over n in Eq. (|14p the Abel-Plana-like summation formula 

y« *ttK) *maf(0) /"".,,, . I™ ,. W) - /(-it) 

This formula is obtained as a special case of the summation formula derived in Ref. [28] by 
using the generalized Abel-Plana formula (see also Ref. |29|). Note that we have the relation 

sin(2A n ) ma 

2A n + (ma) 2 + A2- [U> 

By taking into account Eq. (|17p . we apply the summation formula (|16p with the function 

(18) 



/(z) = y/^ + kj^ + m^a 2 
where we have introduced the notation 



ma 



(ma) 2 + z 2 J ' 



< 9 = < + m z . (19) 
This allows to present the Casimir energy in the decomposed form 



E p+l,q - aE p+l,q + 2E p+l,q + ^ E P+l,q^ ( 20 ) 



where 



4% = -f / (|£r E \A^C> (2D 

is the Casimir energy (per unit volume along the directions z 1 , . . . , z p+1 ) in the topology R p+1 x 
(S* 1 ) 9 when the boundaries are absent. The part 



m N f <ik„ v, / IT /-= T f 00 J^^k\ + m?\ 



n„eZ9 



is the Casimir energy for a single plate (when the other plate is absent) in the half-space. The 
last term in Eq. (|2"pp . 



N f dk„ f°° . \l z2 -K- m l q 



AE p+lq = -- V dz 



n J ( 27r ) P „^. JjK+mg q (z + m)e 2az + z - m 



m 

a(z — m) 



z + m 



, (23) 



is the interaction part. This term is finite for all non-zero distances between the plates and the 
cutoff function can be removed safely. Note that the single plate part of the Casimir energy does 
not depend on the separation between the plates and, hence, will not contribute to the Casimir 
force. 
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The pure topological part (|2"Tj) is investigated in our previous paper [16J. After the renor- 
malization this part is presented in the form 

j=p+2 VD-j+l^j n=1 ' no-^Z ^ 

where we have defined 

D 
l=j+l 

Here and in the discussion below we use the notation 

f v {x)=x v K v {x). (26) 

An alternative expression for the topological part is obtained by using the zeta function technique 
(see below). In particular, the topological part of the Casimir energy is positive for untwisted 
fields [ai = 0) and is negative for twisted fields (a; = 1/2). 



3.1 Single plate part 

Now let us consider the single plate part in the Casimir energy, given by formula (|22p . First 
of all we note that this part vanishes for a massless field. This is directly seen by taking into 
account that in the limit m — > the second term in braces of f|22j) gives nonzero contribution 
which cancels the first term. Another way to see this is the following. For a massless field 
X n = 7r(n — 1/2) and we can apply to the corresponding sum in the Casimir energy (|14p the 
Abel-Plana formula in the form (see, p, 29J) 



E nn - 1/2) = r dx m-i r dx f{ % f ±M , 

n=i Jo Jo e +1 



(27) 



The part of the vacuum energy with the first term on the right of this formula gives the topo- 
logical Casimir energy E^ l q and the second term corresponds to the interaction part AE p+ i tq . 

For the further evaluation of the single plate part for a massive field we apply to the sum 
over n p +2 in Eq. (|22H the Abel-Plana summation formula in the form |30j 

E /(K +2 + a p+2 D = 2 T dx m +iTdx e f }^ + ~J p [T\ • (28) 

n p+2 =-oo - 70 1/0 A=±l 

The part with the first term on the right of this formula gives the Casimir energy for a single 
plate in the case of topology R p+2 x (5 1 ) 9 " 1 and we obtain the following recurrence formula 

c-W - pW _i_ A fW fOQ\ 

E p+l,q ~ E p+2,q-l ' ^P+ZCp+ltf) l Zb V 



where m Ylq _ 1 = yk 2 lq _ 1 + m 2 and we have introduced the vacuum energy per unit volume of 
the compact subspace £ p +i j(J = ^+i q /Vq- I n p9l) - 

A (i) 2NL p+2 ^ cos(27rna p+2 ) ^ [tt 

A p+2 e p+lj9 - -——^-^ („r +2 ) P +2 [ g Wi^mn,.,) 

v 7 y re=l v ^ T y n 9 _ieZ<2-i 



^ m xf p /2 + i(nL p+2 x) 



9-1 



(30) 
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is the part induced by the compactness of the direction z p+2 . In deriving this formula we have 
used the integration formulae 

vrP /2 r((s + 3)/ 2) 



[ dk p r dz(z 2 - k 2 - C 2 )^' 2 f{z) = ^ r (( g + 3 )/ 2 ) r dxx{x 2 _ c 2 )( p+ S +l)/2 /( s 

J P J X /^ K p 1 JKJ T((p + s + 3)/2) J c 1 ; 

(31) 



r (x2 _ 6 2 )(p+ i)/ 2 ^ 2^ 2 + 2 r((p + 3)/2) 
E / b ^ eW^A^.! = E 0F( reLp+2 )P+ 2 cos(2vrna p+2 )/ p/2+1 (nL p+2 6). (32) 

A — zb 1 71 — 1 



and 

( x 2 _ 6 2 )(p +i)/2 - 2f/ 2+2 r((p + 3)/2) 

I ux — — 

/j, p. 

A=±l' 

The first of these formulae is obtained by integrating over the angular part of k p , changing the 

integration variable to y = z 2 — k 2 — c 2 , and introducing polar coordinates in the (|k p |,y)- 

plane. Formula (|32[) is obtained expanding the integrand by using the relation (e u — 1) _1 = 
Y^Li e~ nu and integrating the separate terms in this expansion. 

After the recurring application of formula (|29p the Caimir energy for a single plate is pre- 
sented in the form 

4%, q = V q E% + E^ q , (33) 

where i% is the Casimir energy per unit volume along the directions z 1 , . . . , z D ~ l for a single 
plate in Minkowski spacetime with trivial topology and the second term, 

i=p+2 

is the topological part. The latter is finite and in the corresponding expression the cutoff function 
can be removed. The renormalization is needed for the term Ep ' only. 

3.2 Interaction part 

By using Eq. (|3ip , the interaction part of the Casimir energy is presented in the form 



(4 7 r)-(P+ 1 )/ 2 jV ^ [°° (z 2 -m 2 g )(P +1 )/ 2 
p+1 ' q ~ r((p + 3)/2) n ^ z J mnq Z (z + m)e 2 ^ + z-m 



m 

a{z — m) 



z + m 



(35) 



From here it follows that this part is always negative and it is a monotonically increasing function 
of a. By taking into account the relation 

a(z — m) — m j \z + m) 1 d ( z — m _ 2az 



--In 1 + e~ zaz , (36) 

(z + m)e zaz + z — m 2 dz \ z + m J 

and integrating by parts, Eq. (|35p is written in the equivalent form 

— - e /; - - K i+ ■ <»> 

For a massless fermionic field from here we find 

P+1 ' 9 " - a F( ( p + 3)/2) Al, ^ + 1 



E E^W(«n,), (38) 
n„eZ9n=l 
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where the function f u (x) is defined by Eq. 

Let us consider the asymptotic behavior of the interaction part in the Casimir energy at 
small and large separations between the plates. In the limit L[ 3> a the main contribution 
comes from large values of m, I = p + 2, • • • ,D, and we can replace the summation by the 
integration: ^ n eZ9 — > j dn q . By making use of the integration formula (|3ip with p — > q, we 
find 

A £p+U „ K„ A Ed ,„ = -V, £ ^ - m (l + , (39) 

where AEd,o is the interaction part of the fermionic Casimir energy per unit volume along the 
directions z 1 , . . . , for two parallel plates in D-dimensional space with trivial topology (see 
Refs. [H [23] for the case D = 3 and Ref. [25J for general D). Note that for a massless field we 
have 

A ^,o = - ^^g^ r((g + i)/2) C p + 1), (40) 

where C( x ) * s the Riemann zeta function. 

Now let us consider the limit L\ -C a. In this case and for a; = the main contribution 
comes from the zero mode with n q = and AE p+ x jq /N coincides with the corresponding result 
for the Casimir effect in topologically trivial (p + l)-dimensional space: 

AEp+1 « w iv^ A ^ +1 '° = " r((p+i)/2) ; ' /: : 



xfz 2 



m 2)(p-l)/2 ln A + ^^ e -2aA (41) 

\ z + m J 



where N p = 2^ p+1 - )//2 l. The contribution of the nonzero modes is exponentially suppressed. For 
ai 7^ the zero mode is absent and assuming that am is fixed, to the leading order we have 

A^ +1 ,^- 2(4 ^ )(p+1)/2 c , (42) 

where 

D 

cl=Y. (^fta/Li) 2 , A = min(Q;, 1 - a,). (43) 

l=p+2 

In this case the interaction part of the Casimir energy is exponentially suppressed. 

In the discussion above we have considered the region between the plates. The plates divide 
the background space into three regions: z p+1 < 0, < z p+1 < a, and z p+1 > a. The vacuum 
energy in the regions z p+1 < and z p+1 > a is obtained from the results given above in the 
limit a — > oo. In this limit the interaction part vanishes and we have 

E p+1<q = aEfl lA + E$ ltq , z p+1 < 0, z p+1 > a, (44) 
with the topological and single plate parts given by Eqs. (|24p and (|33j) . 



4 The Casimir force 

The total vacuum energy in the region ^ z l ^ q, I = l,...,p, ^ z p+l ^ a will be 
E p+ i ; qCi ■ ■ ■ c p and the volume of this region is V = c\ ■ ■ ■ c p aV q . The vacuum stress at z p+l = 0+ 
is given by 

P p+ l, q (0+) = -—E p+1 , qCl ■■■c p = P^ q + AP p+hq , (45) 
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where we have introduced the notations 



,(o) _ E p+l,q . „ 15 



?(0) 

^, = AP P+1 , = "^Ai^,,. (46) 

The vacuum stress at z p+1 = a— is given by the same expression. The term fi+i q does not 
depend on the separation between the plates and is the pure topological part of the vacuum 
force. The term AP p+ i^ q is induced by the presence of the second plate and determines the 
interaction force between the plates. Using the formula for AE p+ i^ q , for this part we find 



2(4^)-C+ 1 )/ 2 JV „ f°° zV - mS,)*" 1 " 2 

a^p+i,7 = - r((p+w £ "-g?^-- (47) 



Now we see that AP p +i i9 < independent of the boundary conditions imposed on the field along 
the compactified dimensions and, hence, the interaction forces between the plates are always 
attractive. For a massless fermionic field we have 

AP 2iV ^ ^ ^ n f P /2+i(2ank Qq ) - / p /2+ 2 (2anfc n J 

APp+1 ' q " ^ (2H^ ' (48) 

v ' q n,eZin=l v ; 

For small separations between the plates, L\ 3> a, we replace the summation over n q by the 
integration. In the way similar to that we have used for the Casimir energy, it can be seen that 
in the leading order the interaction force coincides with the corresponding result for parallel 
plates on background of D-dimensional space with trivial topology: 

2N f°° z 2 (z 2 - m 2 ) / 2 - 1 

AP P+1 „ « AP D .„ _ - {4w)D/2rm2) I * W e2 .,' +1 ■ (49) 

The contribution of the nonzero modes is exponentially small. For the massless field we have 

ND(1 - 2~ D ) 

AP ^o = ~ ^^U niD + 1)/2)C(D + 1). (50) 

This result can also be directly obtained from Eq. (|40p . 

For large inter-plate separations, L\ <C a, and for ai = the main contribution comes from 
the zero mode n q = and V q AP p+ i^ q /N coincides with the corresponding result for the Casimir 
effect in (p + l)-dimensional space: 

APp+1 ' q * NpV q APp+1 '° ~ ~ T((p +l)/2)V q L dZ + 1 ■ (51) 

If on ^ and am is fixed the interaction force is exponentially suppressed: 

Nc (p+3)/2 e _2ac 

AP p+1 , q « - {A l a){p+l)/2Vq , (52) 

with Co defined by Eq. (f43l) . 

If the quantum field lives in all regions, in considering the total forces acting on the plate 
we should also take into account the force acting on the sides z p+l = 0— and z p+l = a+. 
The corresponding forces per unit surface are equal to Pjf+i q and they are directed along the 

positive/negative direction of the axis z p+1 in the case P p +i q > 0/P p +i q < 0. Now we see that 
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the topological parts of the force acting from the left and right sides of the plate compensate 
and the resulting force is determined by (|47p and it is attractive for all inter-plate separations. 
There are physical situations [bag model, finite length carbon nanotubes (see below)], where the 
quantum field is confined to the interior of some region and there is no field outside. For the 
problem under consideration, if the quantum field is confined in the region between the plates, 
the total Casimir force acting per unit surface of the plate is determined by Eq. (|45p and the pure 
topological part contributes as well. At large distances this part dominates and the corresponding 
forces tend to increase/decrease the distance between the plates when Pp+i q > ®/Pp+i q < 0- 

In particular, Pp+i q < for untwisted fields and P^+i q > for twisted fields. Hence, if the 
quantum field is confined in the region between the plates, for untwisted fields the Casimir forces 
are attractive for all separations. For twisted fields these forces are attractive for small distances 
and they are repulsive at large distances. 



5 Zeta function approach 

In this section, for the evaluation of the vacuum energy in the region ^ z p+1 ^ a we will use 
the zeta function technique [31] . This allows to obtain alternative representations for the pure 
topological and single plate parts in the Casimir effect. Instead of the divergent expression on 
the right of Eq. (|14p we consider the finite quantity 



E p +i, q (jJ,, s) 



2s+l 



N 



dk p 
(2tt> 



(53) 



n„eZ9 n=l 



where the arbitrary mass scale \x is introduced in order to keep the dimensionality of the ex- 
pression. Performing the integration over k p , we find 



E 



P+1.9 



(54) 



gZ« n=l 



The computation of the Casimir energy requires the analytic continuation of E p+ i tQ (fi, s) to 
the value s = —1/2. The starting point of our consideration is the representation of the partial 
zeta function as a contour integral in the complex plane z: 



1 



d 



(ma I z) sin z + cos z 
1 + ma 



(55) 



where C denotes a closed counterclockwise contour enclosing all zeros A n . We assume that the 
contour C is made of a large semicircle (with radius tending to infinity) centered at the origin and 
placed to its right, plus a straight part overlapping the imaginary axis and avoiding the points 
±OT&n 3 a by small semicircles in the right half-plane. When the radius of the large semicircle 
tends to infinity the corresponding contribution vanishes for Res > (p + l)/2. Assuming that 
(p + l)/2 < Res < p/2 + 1, from (|55p we find the following integral representation for the partial 
zeta function: 



°° -i poo 

E«+ A > 2 ) P/2 " S = -smKs-p/2)]/ dz(z* 
n=l ^ Jm ^ 



d . 
x — ln 

dz 



(m/z) sinh(az) + cosh(az) 
1 + ma 



(56) 
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Hence, the regularized vacuum energy is presented in the form 

(4tt)-p/V s+1 ^ 



2r(s)r(i -s+p/2) 



n„eZ9 " ln « 



d 1 

x — In 

dz 



(m/z) sinh(az) + cosh(az) 



1 + ma 

Now we decompose the logarithmic term in this expression as 



-In 

dz 



(m/z) sinh(az) + cosh(az) 
1 + ma 



d 



d 



a + — ln(l + m/z) H In 1 + 



z — m 



z + m 



-2az 



dz dz 

As a result, we have the following decomposition of the generalized zeta function: 
E P+ i, q {fi, s) = aE^ l q (ii, s) + 2J^2 1>9 (/z, s) + AE p+ i >q ((j,, s), 



(57) 



(58) 



(59) 



where the first, second and third terms on the right hand-side come from the corresponding 
terms in Eq. (|58p . The interaction term AE(fj,,s) in Eq. (|59p is finite at the physical point 
s = —1/2 and gives the result ([3?]) : AE p+ i tq (fi, —1/2) = AE p+ i <q . Below we will be focused on 
the pure topological and single plate parts. 

First let us consider the term (/i, s). This term is the regularized vacuum energy in the 

topology R p+1 x (S 1 )' 1 without boundaries. In this term the integration over z is done explicitly 
and we find 

E (0) , , _ ^NT(s-(p+l)/2) ^ x _ 2s 



(60) 



Further analytic continuation of this expression to the physical point s = —1/2 is done by 
using the extended Chowla-Selberg formula [32J and the corresponding result is given by the 
expression [16J: 



E 



(0) 

P+1.9 



Nm D+1 V q 
(2vr)( s + 1 )/ 2 



W /o J(D+i)/2(mg(L q ,m q )) 



(61) 



m„eZ9 



where we have used the notation 



D 



1/2 



g(L q ,m q ) = \ L 

A=p+2 



2 2 



(62) 



The prime on the summation sign in Eq. (|6ip means that the term m q = should be excluded 
from the sum. 

Now we turn to the part E^^ q (fi, s) which is the regularized vacuum energy in the half-space 
induced by a single plate. In the corresponding integral representation we expand ln(l + m/z) 
in powers of m/z and integrate over z explicitly. This leads to the result 



2s+l N 



8(47r)P /2 r(s) ^ 



, l T((l-p)/2 + s) 



m 



r(z/2 + 



2s-/ 



(63) 



The application to the multiseries over n q of the extended Chowla-Selberg formula allows to 
present E^^ (//, s) as the sum of two parts: 

M /,. „\ _ t/ ^ , fC 1 ^) 



*) = V q E™ (ji, s) + s), 



(64) 
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where 

8(4i)(^r( s )^ ( J r(//2 + i) ' (65) 

is the corresponding quantity in the case of trivial topology R D . The topological term (fi, s) 

is finite at the physical point s = —1/2 and the topological part of the vacuum energy for a 
single plate has the form 

F (M F (i,c), Nm°V ~ 2-n-l) 1 

*p+l, q ~ %+l, q W, W ~ 4(27r) D/2 2L. T (l/ 2 + 1) 

x ^ cos(2vrm 9 • a y )/( / _ D)/2 (m5(L (? ,m (? )), (66) 

m«jeZ9 

where we have used the relation = x~ 2u f v (x). Note that we can write the function 

cos(27rm g -a g ) on the right of formula ([66]) in the form of the product Il£p+2 cos(2-7rmjaj). The 
equivalence of two representations (133p and (j66[) for the topological part in the Casimir energy 
for a single plate can be seen by making use of the relation [16] 

cos(2vrm ? _i ■ OL q ^i)f {l _ D) / 2 {mg{h q ,m q )) 

mq-ieZi- 1 

(2 7 r)fa- 1 )/ 2 L p+2 ^ /(p-0/2+1 (^+2^+2^,,! ) 

W" 1 ^ (m p+2 L p+2 ) P -^2 > ^ > 



n 9 _iSZ9- 



and the formula 



2 f° . c x/ p/2+1 (te) ^ 2-'/ 2 (-l)^ . 

* A - 1 + c 2 = g r(//2 + 1) (bc) ^-0/2 + i(ft). (68) 



valid for sC c < 1. 



6 Special case of topology 

By taking into account the importance of special case p = D — 2, g = 1 in Kaluza-Klein models 
and in carbon nanotubes, in this section we consider it separately. For the later convenience, 
the parameters of the compactified dimension we will denote by = L and «£> = <x The 
corresponding formulae for the separate parts in the Casimir energy take the form 

(o) 2NL-° ^ cos(27rna) 

E d-i,i = (27r) (D+i)/2 2_, n D+i f(D+l)/2(mnL), 



Z)— 1,1 = 2(27r) D/2 cos(27rna) ^ r(f/2 - ^ f(i-D)/2(mnL), (69) 

r((D-i)/2) l f^ oa J mi V ^ + m / 

where we have introduced the notation 

m 2 = [2vr(/ + q) / L} 2 + m 2 . (70) 



iVm D L ^ , ~ 2- / / 2 (-l) i 
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An equivalent representation of the single plate part is obtained from Eq. (|34|) : 



E 



2NL 



(2tt) d / 2 



iE 



cos(2-7rna) 



D 



7T f°° fD/2( n Lmx) 

-f D/2 {nLm) - ^ dx _ = 



(71) 



For the massless case these formulae are simplified to 



E 



(0) 

D-XA 



NL 



AEn- 



D-1,1 



+00 00 

EE 



(2it)- d / 2 N 



n=l 



(-l) n 

2a )D-i ^ ^-^/fl/2(H^ + «IW. 

i=-oon=l 



(72) 



and the single plate part vanishes. In figure [T] we have presented the Casimir energy Ed— 11 for 
a massless fermionic field in the simplest Kaluza-Klein model with D = 4 as a function of the 
inter-plate distance and the length of the internal space measured in units of a fixed length ao- 
The left panel corresponds to the untwisted field (a = 0) and the right one is for the twisted 
field (q = 1/2). For large inter-plate separations the pure topological part dominates and the 
Casimir energy is a linear function of a. At small distances the interaction part is dominant and 
the Casimir energy behaves as a~ D . 




a o E J. 




0.6 

L/a 



0.6 

a/ao 



Figure 1: The Casimir energy of a massless fermionic field in 4-dimensional space with topology 
R 3 x S 1 as a function of the inter-plate distance and the length of the compact dimension. The 
left /right panel corresponds to untwisted/twisted fields. 

In the special case under consideration for the interaction part of the Casimir force we have 
the formula 

2(47T)-( D - 1 )/ 2 iV £5 f°° z 2 (z 2 2 )( D-3)/2 

In the massless case this formula takes the form 

• D ~ 1 ' 1 " "(27TVV2L ^ ^ { } (W 5 \y=4nn\l+a\a/L- (74) 

v ' l=— 00 n=l 
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In figure we have plotted the ratio LAP31/AP3 q versus a/L for different values of the pa- 
rameter a. As it already has been explained before, only in the case of untwisted field the 
Casimir force at large separations tends to the corresponding force (up to the factor related to 
the number of polarizations) for the model where the compactified dimensions are absent. For 
other cases the force is exponentially suppressed at large separations which is clearly seen in 
figure [2j 



3.0 



Figure 2: The ratio of the fermionic Casimir force for two parallel plates in the space with 
topology E? x S 1 to the standard Casimir force in i? 3 , for a massless field, as a function of a/L. 
The values on each of the curves correspond to those of the parameter a. 

7 Applications to finite-length nanotubes 

For a number of planar condensed matter systems the fermionic excitations in the long- wavelength 
regime are described by the Dirac-like model. A well known example is the graphene. In this 
section we specify the general results given above for the electrons on a graphene sheet rolled 
into a cylindrical shape (carbon nanotube). The carbon nanotube is characterized by its chi- 
ral vector Ch = n w &i + m w a2, where ai and SL2 are the basis vectors of the hexagonal lattice 
of graphene and n w , m w are integers. The circumference length of the nanotube is given by 
L = \Ch\ = «g \J n w + m w + n wfnw , with a g = |ai| = |a.2 1 = 2A6A being the lattice constant. 
Zigzag nanotubes correspond to the special case Ch = (n w ,0), and for armchair nanotubes one 
has C/j = (n w ,n w ). All other cases correspond to chiral nanotubes. The electron properties 
of carbon nanotubes can be either metallic or semiconductor-like depending on the manner the 
cylinder is obtained from the graphene sheet. In the case n w — m w = 3q w , q w £ Z, the nan- 
otube will be metallic and in the 7^ 3q w the nanotube will be semiconductor with 
an energy gap inversely proportional to the diameter. In particular, the armchair nanotube is 
metallic and the (n w , 0) zigzag nanotube is metallic if and only if n w is an integer multiple of 3. 

The electronic band structure of a carbon nanotube close to the Dirac points shows a conical 
dispersion EQs) = vp\k.\, where k is the momentum measured relatively to the Dirac points and 
vp represents the Fermi velocity which plays the role of speed of light. The corresponding low- 
energy excitations can be described by a pair of two-component spinors, which are composed of 
the Bloch states residing on the two different sublattices of the honeycomb lattice of the graphene 
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sheet. The corresponding Fermi velocity is given by v F = 3ta/2 (v F « 10 8 cm/s in graphene), 
where t is the nearest neighbor hopping energy. The Dirac-like model is valid provided that the 
cylinder circumference is much larger than the interatomic spacing. For typical nanotubes the 
corresponding ratio can be between 10 and 20 and this approximation is adequate [HE]. In the 
case under consideration D = 2 and we have the spatial topology R 1 x S 1 with the compactified 
dimension of the length L. We will assume that the nanotube has finite length a. As the D = 2 
Dirac field lives on the cylinder surface it is natural to impose bag boundary conditions ([3]) 
on the cylinder edges which insure the zero fermion flux through these edges. The additional 
confinement of the electrons along the tube axis leads to the change of the ground state energy. 
The corresponding expressions for the Casimir energy and force are obtained from the formulae 
of the previous section taking D = 2. Here, by taking into account that in the presence of 
an external magnetic field an effective mass term is generated for the fermionic excitations, we 
consider the general case of massive spinor field. The formulae for a massless case, appropriate 
for carbon nanotubes in the absence of external fields, will be given separately. 

In order to specify the boundary condition on the fermionic field along the compactified 
dimension, W6 not6 that for th.6 {jiwi 

m w ) nanotube the phase factor in the wavefunction has the 
form e l [ m i+( n w-m w )/3]ip ^ w bere ip is the angular coordinate along the compact dimension and mi 
is an integer. From here it follows that for metallic nanotubes the periodic boundary condition 
(a = 0) is realized. For semiconductor nanotubes, depending on the chiral vector, there are 
two classes of inequivalent boundary conditions corresponding to a = 1/3 (n w — m w = 3q w + 2) 
and a = 2/3 (n w — m w = 3q w + 1). In the expressions for the pure topological parts of 
the Casimir energy and force the phase a appears in the form cos(2-7rna) and, hence, these 
quantities are the same for a = 1/3 and a = 2/3. As the boundary induced parts have the 
structure Ylt^-oo f(\l + cx\), the same property holds for these parts. 

In the case D = 2, the general formulae for the separate parts of the Casimir energy from 
the previous section take the form (N = 2) 



1 



—mnL 




n=l 




(75) 



For the Casimir force acting on the edges of the tube we have 




ttL° ^— ' 

n=l 



—mnL 



n 



3 




(76) 
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In the massless case for the total Casimir energy and the stresses we find the formulae 



E 



1.1 



a cos(27rna) 1 



(-iy 



n=l 

oo 



TV 



— y y 

2na ^ ^ n 1 

l=—oo n=l 



fi(Aim\l + a\a/L), 



Pi 



1.1 



1 cos(27rna) 



n=l 



?1" 



2-kcPL ' ~ 2 \y=^n\l+a\a/L- \J7) 



l=—oo n=l 



n- 



The corresponding expressions for the Casimir energy and force in finite length cylindrical nan- 
otubes are obtained from (|77p with additional factor 2 which takes into account the presence 
of two sublattices. In standard units the factor hv-p appears as well. So, for the Casimir force 
acting per unit length of the edge of a carbon nanotube one has: P( CN ) = 2hv F P 1A , where Pi i 
is given by Eq. (|77p . For long tubes, a/L ^> 1, the first term on the right is dominant and 
we have p( CN ) — 0.765hvp/L 3 for metallic nanotubes and p( CN ) m 0.34hvF /L 3 for semicon- 
ducting ones. In the limit a/L <ti 1 the interaction part is dominant. In the leading order the 
Casimir force do not depend on the chirality and one has P( CN > kj -0.lUhv F /a 3 . In fi gure [3] we 
have plotted the Casimir forces acting on the edges of metallic (left panel) and semiconducting- 
type (right panel) carbon nanotube as functions of the tube length for different values of the 
fermion mass. As it is seen, for metallic nanotubes these forces are always attractive, whereas for 
semiconducting- type ones they are attractive for small lengths and repulsive for large lengths. 





Figure 3: The fermionic Casimir forces acting on the edges of the metallic (left panel) and 
semiconducting-type (right panel) nanotubes as functions of the tube length for different values 
of the field mass. 

In the discussion above we have considered bag boundary conditions on the edges of the 
nanotube. The periodicity conditions along the axis correspond to the toroidal compactification 
of the carbon nanotubes. The Casimir energies in toroidal nanotubes are investigated in Ref. 
|16j . where it was shown that the toroidal compactification of a cylindrical nanotube along its 
axis increases the Casimir energy for periodic boundary conditions and decreases the Casimir 
energy for the semiconducting- type compactifications. Recently, in the last paper of Ref. |13j . 
the Casimir interaction between two plates resulting from the quantum fluctuations of the bulk 
electromagnetic field is investigated with one plate being graphene described the Dirac model 
and the other one being ideal conductor. The interaction of the electromagnetic field with the 
fermion field confined on the graphene sheet is equivalent to imposing boundary condition for 
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the electromagnetic field. At large separations the corresponding force is proportional to the 
fine structure constant and falls off as the inverse cube of distance between the plates. 

8 Conclusion 

We have investigated the effect of compact spatial dimensions on the Casimir energy and force 
for a massive fermionic field in the geometry of two parallel plates on which the field obeys MIT 
bag boundary condition. Along the compact dimensions we have assumed periodicity conditions 
([2]) with constant phases ai . The eigenvalues of the wave- vector component normal to the plates 
are roots of transcendental equation (llip . By applying the Abel-Plana-type summation formula 
to the corresponding series in the mode-sum for the vacuum energy in the region between the 
plates, we have explicitly extracted, in a cut-off independent way, the pure topological part and 
the contributions induced by the single plates. The surface divergences in the Casimir energy 
are contained in the single plate components only and the remaining interaction part is finite 
for all nonzero inter-plate distances. The latter is given by Eq. (|37p for a massive field and by 
Eq. (I38p in the massless case. The interaction part of the Casimir energy is always negative. 
We have decomposed the single plate part in the vacuum energy into two terms: the first one 
is the Casimir energy for a single plate in the trivial topology R D and the second one is the 
topological part. The second term is cutoff-independent and in this way the renormalization 
procedure is reduced to that for the plate in topology R D . 

The Casimir forces between the plates have been considered in section HI Single plate parts 
in the Casimir energy do not depend on the plates separation and do not contribute to the 
force. For the region between the plates the forces are presented as the sum of topological and 
interaction parts. In the situations where the quantum field lives on both sides of the plate, the 
topological parts are the same on the left and right sides and the effective force is determined 
by the interaction part only. The latter is given by formulae (|47p and (J48]) for the massive and 
massless fields respectively. With independence of the lengths of compact dimensions and the 
phases in the periodicity conditions, the corresponding force is attractive and is a monotonic 
function of the distance. When the field is confined in the region between the plates only the 
topological part contributes to the resulting force and it dominates at large separations between 
the plates. In dependence of the phases in the periodicity conditions, the corresponding forces 
can be either attractive or repulsive. In particular, for untwisted fields the Casimir forces are 
attractive for all separations and for twisted fields these forces are attractive for small distances 
and repulsive at large distances. For small separations the interaction part dominates and 
the Casimir force is attractive. For small values of the size of the compact subspace and in 
models where the zero mode along the internal space is present, the main contribution to the 
Casimir force comes from this mode and the contributions of the nonzero modes are exponentially 
suppressed. In this limit, to leading order we recover the standard result for the Casimir force 
between two plates in (p + 2)-dimensional Minkowski spacetime. When the zero mode is absent, 
the Casimir forces are exponentially suppressed in the limit of small size of the internal space. 

In section [5] we have evaluated the Casimir energy by using an alternative method based on 
the generalized zeta function technique. With the combination of the extended Chowla-Selberg 
formula, this allowed us to present the topological part for the geometry of a single plate in an 
alternative form given by formula (|66p . As an illustration of the general results, in Sect. [6] we 
have considered a special model with a single compact dimensions. In section [7] we specify the 
general formulae for the model with D = 2. This model may be used for the evaluation of the 
Casimir energy and force within the framework of the Dirac-like theory for the description of 
the electronic states in carbon nanotubes where the role of speed of light is played by the Fermi 
velocity. The pure topological part of the Casimir energy is positive for metallic cylindrical 
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nanotubes and is negative for semiconducting ones. For finite-length carbon nanotubes the 
Casimir forces acting on the tube edges are always attractive for metallic nanotubes, whereas 
for semiconducting-type ones they are attractive for small lengths and repulsive for large lengths. 
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